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Abstract: The determination of a s from hadronic r decays is revisited, with a special em- 
phasis on the question of higher-order perturbative corrections and different possibilities of 
resumming the perturbative series with the renormalisation group: fixed-order (FOPT) vs. 
contour-improved perturbation theory (CIPT). The difference between these approaches 
has evolved into a systematic effect that does not go away as higher orders in the per- 
turbative expansion are added. We attempt to clarify under which circumstances one or 
the other approach provides a better approximation to the true result. To this end, we 
propose to describe the Adler function series by a model that includes the exactly known 
coefficients and theoretical constraints on the large-order behaviour originating from the 
operator product expansion and the renormalisation group. Within this framework we 
find that while CIPT is unable to account for the fully resummed series, FOPT smoothly 
approaches the Borel sum, before the expected divergent behaviour sets in at even higher 
orders. Employing FOPT up to the fifth order to determine a s in the MS scheme, we 
obtain a s (M T ) = 0.320^^, corresponding to a s (M z ) = 0.1 185±g;ooog- Improving this 
result by including yet higher orders from our model yields a s (M T ) = 0.316 ± 0.006, which 
after evolution leads to a s (Mz) = 0.1180 ± 0.0008. Our results are lower than previous 
values obtained from r decays. 
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1. Introduction 

Precision determinations of fundamental parameters within the Standard Model are of 
utmost importance in order to test its internal consistency or point towards physics which 
goes beyond it. In this respect the central parameter of the strong interaction sector is the 
strong coupling a s , and until now tremendous efforts have been put into an ever better 
determination of a s ||, H ■ 

One of the most precise determinations of a s , competitive with the current world 
average, is provided by detailed investigations of the r hadronic width 

„ T\t~ — ► hadrons &v(7)l 

R T = -±- r P^- = 3.640 ±0.010, (1.1) 

F[t -> e v e v T {~j)\ 
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as well as invariant mass distributions |3|, ||, [|, |6], 0]. The recent analyses of the ALEPH 
spectral function data M, ||, |9| on the basis of the final full LEP data set yielded a s (M T ) = 
0.344 ± 0.005 ex p i 0.007th which after evolution to the Z-boson mass scale results in 
a s (Mz) = 0.1212 ± 0.0011. The dominant quantifiable theory uncertainty resides in the 
contribution of as yet uncalculated higher-order perturbative QCD corrections and im- 
provements of the perturbative series through renormalisation group methods. 

Of particular interest for the a s determination is the r decay rate into light u and 
d quarks, R T> v/Ai which proceeds either through a vector or an axialvector current, since 
in this case power corrections are especially suppressed. Theoretically, R T y/A can De 
expressed in the form [||, 10, 11, 12 1 



N c* I t / ,2 



R t ,v/a = -^-^Ew|Kd| 



ud,V/A 
D>2 



r X2) 



Here, Sew = 1-0198 ± 0.0006 @ and 5' EW = 0.0010 ± 0.0010 are electroweak correc- 
tions, comprises the perturbative QCD correction which will be our main interest in this 
work, and the i^L^ denote quark mass and higher D-dimensional operator corrections 
which arise in the framework of the operator product expansion (OPE). The higher-order 
OPE contributions are small and will only be considered towards the end of our work, 
when we present our determination of a s . 

A particular problem emerges from the observation that different ways of performing 
the renormalisation group resummation, namely fixed-order (FOPT) or contour-improved 
perturbation theory (CIPT) JIB], 16], apparently lead to differing results. This is especially 
noteworthy as historically the values of a s extracted from r decays employing CIPT have 
always been on the high side of the world averages, and with the recent update a s (Mz) = 
0.1185 ± 0.0010 [Q] of the latter, this disparity is becoming significant. 

CIPT is conventionally the method of choice, since the expansion of the running cou- 
pling a s (y/s) in a s (M T ) used in FOPT within a certain contour integral (see section |||) 
is near its radius of convergence and thus argued to lead to a poorly behaved fixed-order 
series. This argument, however, is not entirely compelling, since QCD perturbation series 
have zero radius of convergence anyway, and are asymptotic at best, no matter whether 
CIPT or FOPT is used. Indeed, in the large-/?o approximation, which may be viewed as 
a toy model for the entire perturbation series, FOPT was identified to provide the better 
approximation to the full result |jl7|] . With the recent calculation of the 0(af) term in 
the series expansion of the Adler function [18|, the discrepancy between FOPT and CIPT 
appears to be the largest systematic theoretical uncertainty of the a s determination, as it 
is evident that it does not go away by adding the presently known higher-order terms. 

The following study is motivated by the need to resolve this discrepancy and to under- 



stand its origin. Previous investigations jlT], 19 1 show that a preference for CIPT or FOPT 



may strongly depend on the assumptions made on higher-order terms in the perturbation 
expansion. Thus, in section |], we study several toy models in order to address, for each 
model, the following questions: 

i) Are FO and CI perturbation theory seen to be compatible, once terms beyond the 
currently known coefficients of the perturbative series for S^ ' are included? 
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ii) How do FO and CI perturbation theory at a particular order compare to the true 
result for 5^°\ and is the closest approach to the true result related to the minimal 
terms in the respective series? 

hi) And finally, which of the two methods, FOPT or CIPT provides the closer approach 
to the true value at order 0(af ), and in general? 

Here our working assumption is that the true result is approximated with reasonable ac- 
curacy by the Borel sum of the model series, since the power corrections to R T are known 
to be small. 

The lessons learnt from the toy models lead us to proposing an ansatz for the Adler 
function series, more precisely its Borel transform, which incorporates the presently avail- 
able terms in the perturbative expansion, as well as known features of renormalon singu- 
larities p0[| determined solely by the operator product expansion and the renormalisation 
group. This ansatz is described and analysed in detail in section ^. While in the toy models 
discussed in section ||, one may obtain compatible descriptions of the perturbative series by 
FOPT or CIPT, or a preference for one of the two prescriptions, we find that the features 
favouring FOPT prevail in our ansatz for the physical case. For the physical model, CIPT 
never comes close to the result for the Borel sum. On the other hand, FOPT approaches 
this sum in a smooth fashion until its minimal term, after which the expected asymptotic 
(divergent) behaviour sets in. 

We believe that these features are characteristic to R T and therefore argue that FOPT 
provides the better approximation to the perturbative series for both at 0(a*) and 
in general. Based on this observation we proceed to determine the strong coupling a s in 
section |7| in two ways: first, employing FOPT and an estimate of the 0(ot\) term suggested 
by several independent arguments; second, employing our ansatz for the entire series as 
discussed in section ^. Both approaches lead to values of a s (Mz) systematically lower than 
previous determinations from hadronic r decays employing CIPT. 



2. Theoretical framework 



We briefly review the main theoretical expressions required in the analysis of the inclusive 
hadronic r decay width. Further details and complete expressions can be found in the 
original works ||, [2l], 22]. The central quantities in such an analysis are the two-point 



correlation functions 

<^(p) = i I dx<P° <n| t{j%£(z) tfjf (o)t}|n> , (2.i) 



where denotes the physical vacuum and the hadronic vector /axialvector currents are 
given by J^{f(x) = [^7^(75)%] (x). The indices i,j stand for the light quark flavours up, 

down and strange. The correlators n^^(p) have the Lorentz decomposition 

O) = (p^- 9 ,uP 2 )^/ MI \p 2 )+P,P^/ M0 \p 2 ), (2.2) 
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where the superscripts denote the components corresponding to angular momentum J = 1 
(transversal) and J = (longitudinal) in the hadronic rest frame. 

Experimentally, the hadronic decay rate of the r lepton can be separated into the 
contributions of vector R T y and axialvector R Tt A components for the (nd)-quark current 
as well as the contribution with net-strangeness R T ,s> resulting from the (us)-quark current, 



R T = R T ,v + Rt,a + R T ,s ■ 



(2.3) 



In the Cabibbo-suppressed (us) sector, a separation of vector from axialvector contribu- 
tions is problematic, since G-parity is not a good quantum number in modes with strange 
particles. 1 On the theory side, R T can be expressed as an integral of the spectral functions 

^2 



ImnW(s) and Emilys) over the invariant mass s = p of the final state hadrons |23] : 



1 + 2-i_ j Im n« (s) + Im IlW (s) 



(2.4) 



For simplicity, in the following, we shall omit the EW correction factor Sew> but it will 
of course be included in our final numerical analysis for a s . The appropriate combinations 
of the two-point correlation functions resulting from the weak decay through the VF-boson 
are given by 



n( J )(. 



w. 



ud I 



n 



V,(J) 
ud 



(s) + n A J J) 



(s) 



+ \V U . 



(2.5) 



with Vij being the corresponding elements of the Cabibbo-Kobayashi-Maskawa quark- 
mixing matrix. 

The exact (non-perturbative) correlation functions are analytic in the complex s-plane 
cut along the positive axis. Exploiting this property, eq. (^4|) can be expressed as a contour 
integral in the complex s-plane running counter-clockwise around the circle \s\ = || 

2 



R T 



Ml V M r 2 ) 



(2.6) 



|s|=Af2 

The same analytic property holds to any finite order in perturbation theory in a s (although 



the discontinuity is arbitrarily wrong at small s). Eqs. ( |2.4D and (^6) are equivalent if the 
correlation functions are substituted either by the exact values or finite order perturbative 
expansions. The equivalence of eqs. ( |2.4| ) and (2.6) does not hold in renormalisation group 



improved perturbation theory due to the Landau pole singularity [17|. 

Whereas the correlators n^^s) and n ( - * ) (s) themselves are not physical quantities in 
the sense that they contain renormalisation scale and scheme dependent subtraction con- 
stants and thus do not satisfy homogeneous renormalisation group equations, by means of 
partial integration, eq. (|2.6| ) can be rewritten in terms of the physical correlation functions 
D^+°\s) and D^(s), 



D^ 1+0 \s) = 



d r 

Ts 



n ( 1 +°)( f 



£><°>( f 



s d r 
JHds 



: n(°)( f 



(2.7) 



1 A small component with strange quarks but without net-strangeness also resides in R T ,v and Rt,a, 
with the dominant decay channel being t~ — > 7r~ K° K°v t . 
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the first of which being the well-known Adler function 24]. The renormalisation dependent 
contributions to IlW(s) and n^°'(s) drop out after contour integration, and in FOPT the 
perturbative expansions for R T based on II(s) or D(s) are identical. However, after RG 
improvement in CIPT, for the first few terms the series displays a faster rate of convergence 
when employing the correlators D^ 1+0 ^(s) and D^(s). Thus in the present work we shall 
only consider an analysis of R T based on this choice. Utilising the dimensionless integration 
variable x = s/M 2 , eq. (|2.6|) then becomes 



R T 



ITT 



dx 

X 



3 (1 + x) £> {1+0) (M r 2 x) + 4 D (0) (M 2 T x) 



(2.8) 



1*1=1 



For large enough negative s, the contributions to D^ J \s) can be organised in the 
framework of the operator product expansion in a series of local gauge-invariant operators 
of increasing dimension times appropriate inverse powers of s. This expansion is expected 
to be well behaved along the complex contour \s\ = M 2 , except close to the crossing point 
with the positive real axis ||25|| . As can be seen from eq. ( |2.8| ), however, the contribution 
near the physical cut at s = M 2 is strongly suppressed by a zero of order three. Therefore, 
uncertainties associated with the use of the OPE near the time-like axis are expected to 
be very small. Inserting the OPE series for D^ J \s) into (|2.Sj ), performing the contour 
integration, and extracting the terms proportional to |V^| 2 , R T y/A i n the form of eq. ( |L~2^ ) 
emerges. 

The purely perturbative correction 5^ only receives contributions from the vector 
and axialvector correlation function in the chiral limit. Since in this limit vector and axial- 
vector contributions coincide, and D^°\s) = 0, to investigate 5^ we can restrict ourselves 
to the study of the perturbative expansion of the vector correlator Uy (s) in the massless 
case. It exhibits the general structure 

< +0) (*) = - ^ E <L cn, k L\ L = ln-f , (2.9) 

n=0 k=0 " 

with = a(/x 2 ) = a s (fj,)/n and (jl the renormalisation scale. As was already remarked 
above, Uy (s) itself is not a physical quantity. However, the spectral function is as well 
as the Adler function Dy +0 \s), whose general expansion then takes the form: 

n=0 k=l 



In this expression, only the coefficients c n ^\ have to be considered as independent. The 
coefficients c n ^ with k = 2, . . . , n + 1 can be related to the c nj i and /3-function coefficients 
by means of the renormalisation group equation (RGE), while the coefficients c n fi do not 
appear in measurable quantities and c n t n+i 

= for n > 1. Up to order af, the RG 



- 5 - 



constraints lead to: 

C2,2 = - -T Cl,l , C 3i3 = — Ci,i , C 3j 2 = - 7 (P2 Cl,l + 2ft C 2 ,l) , 
C4,4 = - ^ Cl,l , C 4 ,3 = ^ (5/32 Cl,l + 6ft C 2 ,l) , 

C4,2 = - ^ (/% ci,i + 2ft C 2 ,l + 3ft C 3 ,l) . (2.11) 

In our convention, the QCD ftfunction is denned as fta M ) = — {ida^/d/j, = J2k=i PkO>p +l j 
with the first coefficient being ft = HiV c /6 — Nf/S. Since the Adler function Dy +0 \s) 
satisfies a homogeneous RGE, the logarithms in eq. (|2.1C| ) can be summed with the choice 
fj? = —s = Q 2 , leading to the simple expression: 

AT 00 

°r ) w 2 ) = i&£<M"S. (2-12) 



n=0 



where ag = a s (Q)/-K. 

Until recently, the independent coefficients c n ^ were known analytically up to order 



a* At iV c = 3 in the MS-scheme @ they read: 

co,i = ci,i = 1, P2,i = ^ - HCs -[%- fCsW = 1-640, (2.13) 



_ 87029 1103 /- , 275 a ( 7847 262 /- , 25 a \ w , ( 151 19 > \ AT 2 _ n o 71 
C 3 ,l - -nM 4-C3 + -rrC5 - W - TrCs + ttCb -N/ + TM - 27C3 MVf - O- 371 



288 4 1 6 so \ 216 9 ^ 1 9 ^) I 1 \ 162 27 

where numerical values are given at Nf = 3. For the next five- and six-loop coefficients 
C4 t i and C5 i, estimates employing principles of "minimal sensitivity" (PMS) or "fastest 
apparent convergence" (FAC) [^, 30 1, together with known terms of order a^N 2 , exist, 



which for N f = 3 yield (31], |32|: 

c 4 i = 27±16, c 5 ,i = 145 ± 100 . (2.14) 



However, as of this year, the complete result for the 0{oq) coefficient 04,1 is available [18], 
which greatly helps in our analysis. At Nf = 3 it reads: 

_ 78631453 1704247 /- i 4185 ^2 i 34165 a 1995 a _ /in D7fi (11^ 
C 4,l - ^0736 432— « + — C 3 + -gg-C5 ~ — Ql ~ 49.076 . (2.15) 

Since this result turns out to be larger than the estimate presented in eq. (2.14), we shall 



not use the PMS/FAC prediction for 05,1 of ( p. 14 ). Instead, we attempt to estimate this 



coefficient based either on a uniform convergence rate of the series, or on our model. 
3. Renormalisation group summation 

We now discuss the renormalisation group improvement of the purely perturbative cor- 
rection <y(°) to R T by means of resummation of the logarithms appearing in eq. ( 2,1C| ). 
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Returning to eq. (|2.8| ) and inserting the general expansion (2.10) for Dy +0 \s), 5^ is 
found to take the form 



oo n ~\ r A 

* (0) = E<E k ^ h t i (1 - x)3 (1 + x) lnfc " 

n=l fc=i , r , 

\x =1 



-M 2 x 
T 2 " 



(3.1) 



where the identical contribution from the axialvector correlator has already been taken 
into account. 

As discussed above, the Adler function Dy +0 \s) and therefore also 5^ satisfy a 
homogeneous RGE. In fixed-order perturbation theory (FOPT) the logarithms in eq. ( |3.1| ) 
are summed by setting fi 2 = M 2 , leading to 



4o = E^rE^^' 



71=1 



k=l 



The contour integrals Ji in eq. ([T^) are defined by 



Ji 



_ (t —h- x ) 3 (l + x)hi l (-x) = — 



2vu 



2tt 



-fyo + 2 i — 2 I; 3 - J; 4 



(3.2) 



(3.3) 



x =1 



The integrals i/ m are given by 



+- 



/j, m ^(-xr-Mn^-x) 



£ „imo 



-iV +1 

i ( r(/ + 1, —iara) 

m J 



x=l 



(-1) 



l+m 



2l\ 



I(i+l)/2] 



V (-D 



2fc^2A;-l 



m 7T 



fc=l 



(2fc-l)! 



(3.4) 



where r(i + 1, 2) is the incomplete T-function, [n] denotes the integer part of n and m > 1. 
For m = 0, one obtains Ii$ = + (— 1)'] ir l+1 /(l + 1). The first few of the integrals </;, 
which are needed up to order aj, read: 



Jo 



1, Ji 



19 j _ 265 _ 1 „2 t 3355 , 19 2 

12 ' 2 72 3 ' 3 9«R T 19 " 1 



288 1 12 



(3.5) 



in agreement with ref. 16]. 

At order a" FOPT contains unsummed logarithms of order ln*(— x) ~ ir l with Z < n 
related to the contour integrals J/. Contour-improved perturbation theory (CIPT) sums 
these logarithms with the choice /u 2 = —M 2 x in eq. (3.1), which yields 

oo 



(3.6) 



71=1 



in terms of the contour integrals J°(M 2 ) over the running coupling, defined as: 

J a n (M 2 ) = JL j ± (l _ x f (l + x ) a n (-M 2 T x) . (3.7) 



1*1=1 
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Figure 1: Contour integrals J^(M^) of eq. (3.7) required for CIPT as a function of the perturbative 
order n, computed with an initial value a s (M T ) = 0.34 and 4-loop running. 



In contrast to FOPT, for CIPT each order n just depends on the corresponding coefficient 
c n> i. Thus, all contributions proportional to the coefficient c n> i which in FOPT appear at 
all perturbative orders equal or greater than n are resummed into a single term. This is 
related to the fact that CIPT resums the running of the QCD coupling along the integration 
contour in the complex s-plane as can be derived directly from eq. ( |2.12 ). 



In view of our numerical analysis below, a few additional remarks are in order. Al- 
though from the form of eq. (|3.6|), the contour integrals J^(M^) could be considered as 
effective couplings, they have a rather non-trivial dependence on the perturbative order n. 
This can be seen from figure [j], m which we display J^(M^) / a(M^) n for an initial value 
a s (M T ) = 0.34 Q as a function of n up to n = 30. Up to the 7th order, J^(M^) is positive 
and then above the 7th order turns negative. This implies that Jj(M^) is rather small 
which is later reflected in the fact that with four-loop running of a s always the 7th term 
in the CIPT series is found smallest. This observation already casts some doubts on the 
approach of treating the CIPT series in the sense of an asymptotic series for which quite 



often the optimal truncation is provided by breaking the series at the smallest term [33|. 

We now recall (as is well-known, see for instance (]l8| ) that the two approaches lead to 
significant numerical differences. Using the analytically known coefficients of eqs. ( 2.13; ) as 



well as (2.15) and a s (M T ) = 0.34 in eqs. (|3.2|) and (p^q), we obtain: 



a] a 2 s a 3 s aj a 5 s 



5^ = 0.1082 + 0.0609 + 0.0334 + 0.0174 (+ 0.0088 ) = 0.2200 (0.2288) , (3.8) 
<@ = 0.1479 + 0.0297 + 0.0122 + 0.0086 (+ 0.0038 ) = 0.1984 (0.2021) . (3.9) 



The CI series displays a faster convergence, but the two series do not appear to approach 
a common value as successive terms are added. Summing both series up to order a^, the 



difference between FO and CI perturbation theory amounts to 0.0216. The size of this 
difference is of the order of the last included term in the FO series and about a factor 
of 2.5 times the corresponding CIPT term. This apparent disparity in the perturbative 
prediction at the moment represents the dominant theoretical uncertainty in the extraction 
of the strong coupling a s from the hadronic r decay rate. Thus, it is legitimate to ask 
how the series for <5qj and 5^q behave if even higher-order perturbative coefficients are 
included. On the one hand both results would be expected to be compatible, if an all-order 
result were available. On the other hand, re-expansion of the contour-improved integrals 
J"(M 2 ) into the fixed-order series in a(M^) results in a series that is barely convergent for 
realistic values of a(M 2 ) p6| , so one might question the validity of FOPT altogether. In the 
following sections we examine these issues assuming different behaviours of the higher-order 
terms in the series. 

An immediate question which arises is, if, on the basis of the series up to the fourth 
order, we are in a position to say something about the next coefficient 05,1? As a somewhat 
naive guess, we might assume that the size of the fifth order should be at most of the size 
of the previous term and larger than zero, which is to say that the asymptotic (divergent, 
sign-alternating) behaviour has not yet set in at the fifth order. Applying this criterion to 
the CI series, one arrives at the estimate < c$ ; i < 642. A slightly more elaborate estimate 
can be based on the striking feature that the convergence rate of the FOPT series is found 
to be very uniform. Each new term is rather closely half of the preceeding one, and the 
slight dependence on the order can even be nicely fit to a linear behaviour. Assuming that 
this property persists also at the fifth order, we arrive at the estimate 2 

c 5 ,i « 283. (3.10) 

This value lies close to the centre of the range given above, and in section |6] it will be seen 
that eq. ( |3.10|) is corroborated by our model of higher order coefficients. Interestingly, it is 
also close to the update of the FAC estimate to account for the newly available exact C4 1, 
which yields c 5j i = 275 @. Including the estimate (|3.10| ) in the series for 5^q and , 
the numbers in brackets given in eqs. (|3.8| ) and (|3.9| ) are obtained. Now, the difference 
<5p°Q — 6qi = 0.0267 is increased even further, and found much larger than the last included 
summands. 

4. Higher orders: toy models 

To acquire some feeling of what can happen to 5^ in FO and CI perturbation theory when 
higher terms in the perturbative series are included, in this section we exhibit a few toy 
models. Let us emphasise that we do not believe that these models have much in common 
with the true QCD case (with the exception, perhaps, of the large- 0q approximation). 
Rather our concern is to find out which features of the higher-order series determine whether 
CIPT of FOPT represents a better approximation to the true result. Inspired by what we 

Postdicting 04,1 in this way, we find 04,1 = 52 ! The first few terms of the FO series for tf<°> are very 
nearly geometric. 
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Figure 2: Results for (full circles) and (grey circles) at a s {M T ) 
order coefficients c n _\ = for n > 6. The results up to the fifth order coincide with eqs. (3 
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Perturbative order n 

(0) 



18 



0.34, setting the higher- 
and 



(3.9) respectively. 



learn from the models, we are led to the construction of a realistic ansatz which will be 
discussed in detail in the following sections. 

4.1 Truncated Adler function 

Let us begin by considering the case in which all perturbative coefficients c U) i are set to 



zero for n > 6. This model has already been investigated in ref. [19|, in order to see 
how FOPT and CIPT might be compatible even though at the fifth order their difference 
appears rather dramatic. A graphical display of the model is shown in figure ^. The result 



for dpi is given as the full circles and as the grey circles, as a function of the order 



t(0) 



^FO 



CI 



up to which the perturbative series has been summed. To guide the eye, we have also 
connected the points by straight line segments. Since at each order 5q) only depends on 
the explicit coefficient c^i, starting from the order where the coefficients are taken to be 
zero, CIPT becomes exact (and thus 5qj constant). On the contrary, in FOPT, besides the 
contribution from the c nj i, at each order we also have contributions involving all lower i 
with k < n, which are due to the running of a s along the complex contour. In FOPT these 
latter terms are present even if the higher c n) i are set to zero, and entail that for higher 



orders S^q oscillates around the constant Sq^ . It is evident that for the present example 
FOPT represents a rather poor approximation to the exact result up to very high n. 
Quite generally, we can write <5pQ in the form 



t(o) 



5 (o) 

J FO 



2\n 



(4.1) 



n=l 



where the c n ^\ series is simply the Adler function series, while the g n series represents 
the additional contribution from the contour-integral of the Adler function series. By 
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comparison with eq. (fT^), g n = Ylk=2 k c n,kJk-ii hence depending only on c^i with k < n 
and /3-function coefficients as remarked above. With only a few values of c n ^ given, the 
g n series has a finite radius of convergence [16]. For a s (M T ) slightly larger than 0.34 the 
series becomes divergent due to large running coupling effects along the circular contour - 
the amplitude of oscillations of FOPT around the exact result grows and FOPT is never a 
good approximation. This is the reason why CIPT is usually argued to provide the more 
reliable approximation to R T . 

The general picture observed in this special case also translates to other models in 
which the coefficients c Ut i are small compared to the running effects along the complex 
contour. However, the real series expansion of R T is not of this general form. Rather, both 
the c n and g n series are divergent with zero radius of convergence, and can at best assumed 
to be asymptotic. Moreover, systematic cancellations are predicted between the c n and g n 
terms, when n is sufficiently large. 

4.2 The "large-/3o" approximation 

We take a first look at the issues that arise when the series have zero radius of convergence 
in another toy-model, the so-called "large-/?o" approximation. 3 An analytic result for the 



Borel transform of the Adler function and the corresponding R T is available |17j, f2Q, 34 



35, 36] in this approximation, and thus the perturbative coefficients c nj i are known to all 
orders. 

Let us briefly review the results for 5^ in the large-/3o approximation. To make contact 
to the notation employed in the original works on renormalons in connection to the Adler 



function and R T [17, g(| (for a review see |p0[), it is convenient to define the new function 
DO) by 



U7 L D $ +0 \s) 1 + D(s) ee l + J2r n a s (V~s) n+1 . (4.2) 
,(i+o). 



n=0 



Then the expansion coefficients of D v (s) and D(s) are related by c n> i = ir n r n _i. Next, 
the Borel-transform of D(s) is defined by 



N t 



B[D](t) ee £r„-. (4.3) 



n! 

n=0 



If B[D](t) has no singularities for real positive t (which is not the case for the Adler 
function) and does not increase too rapidly at positive infinity, one can define the Borel 
integral (a positive) as 

oo 

D(a) ee J dte-V a B[D](t), (4.4) 
o 



3 For historical reasons, we shall speak about the "large-/3o" approximation, although in the notation 
employed in this work, the leading coefficient of the /3-function is termed f)\. The "large-/3o" approximation 
uses only the term with the highest power in the number of light flavours, Nf, in c„,i, and replaces Nf by 
— 3/3i . Correspondingly, in the evolution of a s , only one-loop running is taken into account. 
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which has the same series expansion in a as D(s) does in a s (y/s). The integral D(a), if 
it exists, gives the Borel sum of the original divergent series. Calculating so-called bubble- 
chain diagrams, it was found that the Borel-transformed Adler function B[D](t) obtains 
infrared (IR) and ultraviolet (UV) renormalon poles at positive and negative integer values 
of the variable u = /3\t /(2tv), respectively [34, |35|. (With the exception of u = 1.) While 
the IR renormalons are related to power corrections in the operator product expansion, 
the leading UV renormalon, being closest to u = 0, dictates the large-order behaviour of 
the perturbative expansion. 

The central result of refs. |34|, 35] is that the Borel transform of the Adler function in 



the large-/3o approximation (see also eq. (5.10) of ref. |2(|) can be expressed as [35] 

_ p-Cu °° (—1\ k k 

(4 - 5) 

where C is a scheme-dependent constant which cancels the scheme dependence of a in 
eq. Q4.4D , such that D(s) is independent of this choice. In the MS-scheme it takes the 
value C = —5/3. Taylor expanding the Borel transform in the variable u and performing 
the Borel integral (<L4) term by term, the perturbative coefficients c n .\ for the large- (3q 
approximation can be deduced. Numerical values for the first 12 coefficients in the MS- 
scheme are presented in table |]. One observes that the dominance of the leading UV 
renormalon at u = —1, and the corresponding sign alternating behaviour of the series, in 
the large-/?o approximation sets in at the 6th order. 

The perturbative coefficients of table [l] can now be employed to calculate 5^ in the 
large-/?o approximation. For consistency the running of a s was also implemented at the 
one-loop order. (We shall soon see that this plays an important role.) A graphical account 
of our results at the physical coupling a s {M T ) = 0.34 is displayed in figure ||[ Again, the 
full circles correspond to 5^q in FOPT while the grey circles provide 5q) in the large- 
st) approximation. In both cases, the grey diamonds represent the order at which the 
FO and CI series have their smallest terms, before the asymptotic behaviour sets in. The 
qualitative picture is rather different from the previous case of the truncated Adler function 
displayed in figure ^. Both series appear to reach a plateau before the divergence sets in, 
though the plateau is reached at higher order for FOPT in agreement with the earlier 
analysis [17]. However, the difference in the plateau values of 5^ is far larger than the 
minimal terms of both series, hence FOPT and CIPT seem to give incompatible results 
within the conventional uncertainty estimates. 

4 A particularly effective way of analytically generating the large-/3o coefficients c n ,i can be derived from 
eq. (26) of ref. |7|. 



ci,i 


C2,l 


C3,l 


C4,l 


C5,l 


03,1 


1 


1.5565 


15.711 


24.832 


787.83 


-1991.4 


C7,l 


C8,l 


C9,l 


ClO,l 


eii,i 


Cl2,l 


9.857 • 10 4 


- 1.078 • 10 6 


2.775 • 10 7 


- 5.388 • 10 8 


1.396 • 10 10 


-3.598 • 10 11 



Table 1: Perturbative coefficients c ni i in the large-/?o approximation up to 12th order. 



- 12 - 



0.25 



o 

"to 




0.15 



Figure 3: Results for <5pQ (full circles) and 5^ (grey circles) at a s {M T ) — 0.34, employing the 
higher-order coefficients c n ^\ of table [l] obtained in the large-/?o approximation, as a function of 
the order n up to which the terms in the perturbative series have been summed. The straight line 
represents the result for the Borel sum of the series, and the shaded band provides an error estimate 
inferred from the complex ambiguity. 



6 8 10 
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Which is the better approximation to the "true" result? Since we now have a closed 
expression for the Borel transform of the Adler function, we can compute the Borel integral 
and compare the perturbative series to this result. We do not expect the Borel integral 
to correspond exactly to the "true" result. First, the "true" result receives condensate 
corrections from higher-dimensional operators in the OPE expansion. Second, when the 
Borel transform has poles on the positive i-axis, the Borel integral must be defined by an 
arbitrary deformation of the contour into the complex plane, which introduces an ambi- 
guity, whose dependence on a matches the condensate corrections. The two are closely 



related (see the review [p0| ), and it has been observed [17| that the size of the ambiguity 
divided by ir is indeed of the order of non-perturbative corrections in the OPE. We thus 
conclude that we expect the "true" result to coincide with the principal value of the Borel 
integral within an accuracy set by about the ambiguity of Borel integral (divided by ir), 
and certainly not parametrically larger. 

In figure ^ the horizontal line represents the value of the Borel integral when the series 
is summed via the principal-value prescription, and the shaded region marks the size of 
the complex ambiguity in the Borel integral. 5 Details on the analytical calculation of the 
Borel transform have been relegated to appendix A. A good approximation of the "true" 
result should approach the shaded region smoothly and diverge eventually. Figure || shows 



the remarkable result, already observed in [17], that in the large-/?o approximation FOPT 
approaches the Borel sum in a rather monotonous fashion until the 10th order after which 
the sign-alternating divergent behaviour of the series sets in, while CIPT never comes close 



5 To obtain the maximal complex ambiguity the IR renormalon poles are circled in such a way that they 
all contribute with the same sign. The modulus of the imaginary part thus obtained is then divided by n. 



- 13 - 



to the "true" result due to an earlier onset of the divergence. Clearly FOPT represents the 
better approximation here, even at n = 4, in stark contrast to the previous example of the 
truncated Adler function series. 

Let us try to gain a better insight as to why the behaviour in the large- (3q approximation 
is so different from our first model. An immediate observation that can be made is that the 
c n> i have to be of a similar size as the running effects g n from the contour, such that strong 
cancellations between these two effects take place. To make this more precise, consider 
the simple relation between B[D] and the correspondingly defined Borel transform of 
expanded in a s (M T ) given by |p8| 



B[5 {0) ]{u) = B[D](u) sin(7ru) 



1 



+ 



f 



TTU 7r(l — u) 7r(3 — u) 7r(4 — U) 



(4.6) 



which is valid in the large-/3o limit. 6 We are then in a position to derive the large-order 
behaviour. Applying the decomposition (|4.1|), we obtain 



9n+i = 1 y J n\ 



Pi 



-e" 5/3 (-l) n 
9 V ; 







e 10/3 


n + 


D 


H 




2 n 






e 10/3 


n + 








2 n 



+ ... 



+ ... 



(4.7) 



where the dots denote the less important poles beyond the leading ultraviolet renormalon 
pole (u = —1) and the leading infrared one (it = 2). This shows explicitly the strong can- 
cellations that take place between the Adler function series and the extra terms generated 
by the integration along the circle. The suppression of the large-n divergence allows FOPT 
to approach the Borel sum smoothly. We also see why CIPT fails in this case: when such 
strong cancellations between the c n ^\ and g n series are present, it is mandatory to combine 
the two series at the same order. However, CIPT uses the c n ^\ up to some finite order, 
while summing the g n to all orders, thus missing the cancellation [17, pOf |. The result is 
that CIPT runs earlier into the leading UV renormalon divergence, as seen in figure |3|, 
though the divergence is damped by the suppression of the effective couplings J^(M^) as 
discussed in section ||. 

There are some lessons that can be drawn from our observations which are valid beyond 
the large-/?o toy model. First, since the known exact coefficients of the Adler function series 
show no interference of a sign-alternating component, we expect the leading IR renormalon 
at u = 2 to be the most relevant contribution, before the eventual sign-alternation takes 
over, just as above. Second, the leading IR renormalon contribution will no longer cancel 
completely as it does in eq. ([O]). However, it is true in general that it is suppressed by 
a factor 1/n 2 in the sum c ni \ + g n relative to c n> i alone pq ], since this follows from the 
OPE and the anomalous dimension of the gluon condensate. Thus, we expect some of the 
features of the large-/?o toy model to survive in the realistic QCD case. 



6 To derive this result, insert eq. (4.4) into eq. (2.8) and perform the contour integral employing the 
one- loop expression for Q s (y / s). 
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Next, we corroborate these findings with two further simplistic models, which will also 
uncover another crucial ingredient for the building of a more realistic model for higher 
order coefficients. 



4.3 Single pole models 

In order to separate the effects of a specific renormalon pole, and to give support to our 
previous findings, we investigate a single IR renormalon pole at position p, for which we 
shall consider the particular cases p = 2 and p = 3. Explicitly, our single pole model for 
the Borel transform of the Adler function series takes the form: 

d m 3 

B[D p ](u) = _£_ + ^dfO„*. (4.8) 
VP u ) i=0 

We fix the residue such that the perturbative coefficient = 283, and the polynomial 
terms df are adjusted such as to reproduce the lower order coefficients according to 
eqs. flOj) and ( grp . The reason behind employing C5 1 to fix the residue is that we like 



to work with a perturbative order at which a dominance of the leading IR renormalon to 
the coefficient is expected, but this assumption is not crucial for our argument. 

The numerical results for the model ( [4.8D with p = 2 are shown in figure ||. The upper 
plot corresponds to SpQ, the lower one to 5^. Let us first concentrate on the triangles 
which are connected by dashed line segments, and for which the running of a s at 1-loop 
has been employed. This case is rather similar to the large-/3o approximation. FOPT 
approaches the Borel sum (the horizontal dashed line) well, although here in a slightly 
oscillatory manner, and up to the order shown in the plot the divergent behaviour has not 
yet set in. Conversely, like for the large-/3o approximation above, for CIPT the divergent 
behaviour sets in much earlier and the series never comes close to the Borel sum. 

The picture drastically changes when the running of the coupling in the contour in- 
tegration is implemented at four loops. This case is shown as the full and grey diamonds 
in figure which are connected by dash-dotted lines. The corresponding Borel sum is the 
straight dash-dotted line. 5^q first overshoots the "true" result by a large amount, then dis- 
plays a large oscillation similar to figure ||, before it starts to diverge, while CIPT appears 
more like the 1-loop case. Still, with 4-loop running, neither FOPT nor CIPT approach the 
Borel sum in a sensible fashion. The reason for this unexpected behaviour can be traced 



back to the fact that our model (4.8) only contains a simple pole. Such a simple pole is the 
correct structure of a renormalon pole in the large-/?o limit, but when higher terms in the 
/3-function are to be included, the renormalon pole structure gets more complicated, also 
involving cuts, whose gross features are determined by the OPE pC| ]. Hence, to construct 
a consistent model which aims to use a 4-loop running coupling, also the renormalon cut 
structure has to be incorporated at the same order. The results required at four loops 
will be derived in the next section. Once four-loop running is consistently included in the 
Adler function Borel transform eq. (fl^) and the contour integration, we obtain the circles 
connected by solid lines in figure |4j. Now we find again that FOPT for an IR pole at u = 2 
smoothly approaches the Borel sum, while CIPT fails. 
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Figure 4: Results for S^q (full circles, diamonds, triangles) and <$£?j (grey circles, diamonds, 
triangles) at a s (M T ) = 0.34, for the model of eq. (4^) at p = 2. The triangles correspond to a 
running of a s at 1-loop, while for the diamonds a 4-loop running coupling has been employed. For 
the circles the pole in eq. (4.8) is modified to account for the 4-loop renormalon cut structure. The 



straight dashed, dash-dotted and solid lines represent the respective Borel sums. 



To conclude the discussion of simple toy models, we finally investigate the ansatz fld.gj) 
with p = 3. Graphically, this case is displayed in figure ||, with the same notation as used in 
figure We see that this case very much resembles the model of figure where the higher 
c n> i had been set to zero. The p = 3 model differs from the p = 2 one in two respects. 
First, it follows from eq. ( |4.6j ) that there is no suppression of c n .i + g n relative to c n) i in 
large orders. More importantly, with p = 3 the divergence of the series is milder. Hence, 
the Adler function coefficients c n ^ are much smaller than the running effects g n along the 
complex contour, an expectation that can be verified by explicitly investigating the c nj i in 
this model. Thus, like for the first model, the truncated Adler function, CIPT provides a 



- 16 - 



0.24 



0.22 



%° 0.2 



0.18 




Borel resummation 
"Hoop running, 1 -loop pole 
4-loop running, 1-loop pole 
4-loop running, 4-loop pole 





10 12 
Perturbative order n 

Figure 5: Results for SpQ (full circles, diamonds, triangles) and S^i (grey circles, diamonds, 
triangles) at a s {M T ) — 0.34, for the model of eq. (4^?) at p = 3. The triangles correspond to a 
running of a s at 1-loop, while for the diamonds a 4-loop running coupling has been employed. For 



the circles the pole in eq. (h8) is modified to account for the 4-loop renormalon cut structure. The 
straight dashed, dash-dotted and solid lines represent the respective Borel sums. 



good account of the Borel sum, and FOPT is only able to approach this value with large 
oscillations. 

4.4 Resume 

The main conclusions from the toy examples are as follows. We find that CIPT provides 
the better approximation whenever running coupling effects to the series expansion of R T 
dominate over the intrinsic Adler function coefficients as should have been expected. This 
is the case in truncated perturbation theory and in models with weak factorial divergence, 
such as the p = 3 single-pole model. We find that FOPT provides the better approxima- 
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tion, whenever there are systematic cancellations between the Adler function and running 
coupling contributions. Such cancellations occur in the p = 2 single-pole model, the large- 
/3q approximation, and in general for the leading IR renormalon contribution to R T . We 
also find that to correctly account for these cancellations, the running coupling effects 
have to be implemented at the same loop-order in the contour integral and the renormalon 
structure of the Adler function Borel transform. The interesting question is now which of 
these features is relevant to the "real world". 

The discussion of this section already allows to disfavour a large class of models for the 
realistic Borel-transformed Adler function that were initially considered by us. In this class 
fall models with renormalon poles that only have integer power, because they fail to account 
consistently for running coupling effects beyond 1-loop evolution. For the same reason - 
besides not using the available information on the known positions of the renormalon poles 
- also models based on Pade approximation are of limited use. 

In order to be able to build a more realistic model for B[D](u), as a prerequisite, 
employing the RGE and the structure of the OPE, in the next section we shall derive the 
general form of the renormalon cut including /^-function effects up to 4-loop level. 



5. Renormalon poles at four loops 

Our next aim is a physically motivated ansatz for the Borel transform of the Adler function 
B[D](t), which incorporates all known exact results, and on the basis of which we will be 
in a position to investigate the influence of higher-order perturbative contributions, and to 
perform a comparison of the perturbative series in FOPT and CIPT with its Borel sum. 
As has been seen in the last section, if higher-order running effects are to be included in 
the contour integration, the renormalon pole structure should match the corresponding 
loop order. The derivation of the renormalon cut is based just on the structure of the 
OPE and the RGE, detached from any limitations of the large-/3o approximation, up to 
an unknown overall constant. The expressions obtained in this section extend the analysis 



already presented in sections 3.2.3 and 3.3.1 of ref. [20] to one more loop order. 

Let us begin with the IR renormalon poles. The central idea is that the IR renormalon 
ambiguity of the Borel integral arises from long-distance regions in Feynman integrals, and 
therefore must be consistent with the power-suppressed terms appearing in the operator 
product expansion [39, jicR . Comparing the energy dependence of a certain term in the 



OPE to the one of the complex ambiguity of the Borel integral, the renormalon singularity 
that gives rise to this ambiguity can be determined. A generic term in the OPE of D(s) 
from an operator Od of dimension d can be written as 



c O d ( a Q) 



(Pa 



[°q]' 



J*) 

0i 



/=J(0) , fr(l) „ , A(2) 2 , 
L O d + L O d a Q + L O d a Q + 



Q a 



(5.1) 



where the anomalous dimension 70 d of the operator Od is defined by 



J 1 ) a +J 2 ) a 2, (3) 3, 



O d {ji) • (5-2) 
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For convenience, we have expressed eq. ( |5.l[) in terms of the scale invariant operator Od, 
defined by 

" lo d { a ») 



O d = Odifx) exp 



P> M ) da » ( ' 



(5.3) 



such that higher order coefficients of 70 d are contained in the Wilson coefficients Cq . 
Since eq. (5J.) will only be needed up to a multiplicative factor, we do not have to specify 
the constant of integration in (|5.3|) , and without loss of generality, it can be assumed to be 
zero. Employing the RGE for qq, the Q-dependent part of ( |5.1[) can be written as: 



c oA a Q) 
Q d 



const, x C 0d (a Q )e ^ a Q [a Q ] i expi d 



1 + h 



(3{a) (3 ia 2 (3fa 



da 



-d 



#2 r 



const, x C 0d (a Q ) e ^ a Q [a Q ] ~ B 'i 1 + b 1 a Q + b 2 a 2 Q + . 



(5.4) 



where the coefficients b\ and b 2 are found to be: 

1,2 



ft i = i ~ Pifo) » 6 2 = 77 - t4 (p1 ~ 2 frfo & + / W • (5.5) 
Pi ^ ^Pi 



To find the Borel transform that matches the Q-dependence of (5^4), we take the 
ansatz: 



B[D™](u 



dp 



1 + h(p - u) +b 2 (p - uf + 



(5.6) 



(p-u) 1+ "f 

Employing eq. ( |A,8j ), the imaginary ambiguity corresponding to the Borel integral of 
B[D p K ](u) is found to be: 



Im 



_2p_ 



D Lii (a Q ) = const, x e 01 a Q [c 



1 + h — 7 a Q + b 2 ^- 7(7 - 1) a 2 Q + . . . 



(5.7) 



Comparing eqs. fl5.4|) and (|5.7D, one deduces: 



d 



,(l) 



P = t: > 7 = 2p-9 - —f- , bi 



2(6i +ci) ~ 4(6 2 + 6ici + c 2 ) 
, »2 = 



ft 2 7(7-1) 



2 ' # pi ' P17 

< JC 0d and c 2 = C 0d /G 0d 
performing the Borel integral term by term yields the perturbative series: 

tt<£ r 

4- 

n=0 



(5.8) 



where c\ = Cq'/Cq' and c 2 = Cq 1 /Cq 1 . Taylor expanding the ansatz fl5,6| ) in u and 



pIR( a ) = J d f y r(n + 1 + 7) ( ^Xa n +1 



x 2p (6 1 +ei) + /2pV (6 2 + ^iCi + c 2 ) + ^ / 1 



/3i (n + 7) V/V (n + 7)(n + 7-l) 



(5.9) 



Eq. ( |5.9| ) extends the corresponding eq. (3.51) of ref. |20| to include terms of order 1/n 2 in 
the large-order behaviour of the perturbative series. 



7 Note a missing sign in the global factor containing /3o in eq. (3.51) of 120] , which should read (— 2[3o/d) r ' 
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The corresponding expression for a general UV renormalon pole when higher orders 
in the running are included can be deduced by formally considering negative coupling oq 
and a Borel integral that ranges from zero to minus infinity flil] ]. Then one finds poles 
leading to ambiguities on the negative real axis, and the complex ambiguities can again 
be identified with the RGE properties of some higher-dimensional operators. The result is 
that the constants in the parametrisation of the general UV singularity, 



B[D^](u) 



(p + u) 1+ 7 



1 + h(p + u) + b 2 (p + uf 



(5.10) 



can be obtained from the corresponding parameters of the IR renormalon pole ( |5.6| ) with 
the replacement p — ► — p, leading to: 



- - tt + ~ + • h = - bl\p ->■ -p] , b 2 = b 2 [p -> -p] ■ 
The perturbative expansion corresponding to eq. ( 5.10| ) takes the form: 8 



(5.11) 



pi+7T(l + 7) 
1 + &1 



Er(n + H7)(-|]C 



n=0 

P7 



(n + 7) 



+ b 2 



(n + 7) (n + 7 — 1) 



(5.12) 



In the next section, the general forms of the IR and UV renormalon singularities 
(|5.6| ) and ( [5,10| ) will be employed to construct a physically motivated model for the Borel 
transform of the Adler function B[D\(u). In order to describe the leading IR renormalon 
at u = 2 as well as possible, in this case we include the known Wilson coefficient function 
and anomalous dimension of the gluon condensate (aG 2 ) ^2|, which results in [43] 



Ca _ Cf_ 
2 4 



fh. 
A' 



Cl [(aG 2 )] 

while all other c's appearing in the equations above will be set to zero. 



(5.13) 



6. A physical model for the Adler function series 

To clarify whether FOPT or CIPT results in a better approximation to the r hadronic 
width, we need to construct a physically motivated model for the Adler function series 
beyond the order n = 4, up to which it is known exactly. Since the four lowest coefficients 
are available, it is reasonable to attempt to merge the low-order series with the expected 
large-order behaviour. Furthermore, it is convenient to generate the series from its Borel 
transform. 9 A physical model should account for the following features: 

8 There is a sign mistake in the term proportional to 1/n in the corresponding eq. (3.48) of J2(J. 
9 A previous attempt to resum the perturbative series for R T , based on Borel transforms of the Adler 
function, has been made in ref. J44]. 
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It should reproduce the exactly known c nj i, n < 4. 

The very large order behaviour is governed by a sign-alternating UV renormalon 
divergence. Since the low-order series shows no sign of a sizeable alternating com- 
ponent, it should be sufficient to include the leading singularity at u = — 1. On the 
other hand, since the intermediate orders are governed by IR renormalons, at least 
two IR renormalon singularities, at u = 2 and u = 3, should be included to merge 
the large-order behaviour with the low-order exact coefficients. 

Since four-loop running is employed in the contour integral that relates the Adler 
function to R T , the renormalon singularities in the Borel transform cannot be simple 
poles, but should be extended consistently to cuts. This is particularly important 
for u = 2. We use the equations from section [|, incorporating c\[{aG 2 )\ into the 
description of the u = 2 cut. For the UV renormalon singularity at u = — 1, we use 



7 = 1 — 2(32 /P( in eq. ( 5.11 ), since it is known to be a double pole in the large-/3o 



approximation. 10 
We are thus led to the ansatz: 

B[D](u) = B0Y y ](u) + B[D l 2 R ](u) + B[Df](u) + d%° + df°u, (6.1) 
where the first three terms use ( |5.6[ ) and ( 5.10| ) as building blocks for the three leading 



renormalon singularities. 11 The model (O) depends on five parameters, the three residua 
of the renormalon poles, d^ v , and d^, as well as the two polynomial parameters d^ 
and di° . To fix these parameters, we match the perturbative expansion of our model to 



the known coefficients c nj i of eqs. pl3|) , ( |2T5| ) and ( pUPp . First, the residua are fixed 



using the coefficients 03,1, c^i and 05,1, and then the polynomial coefficients are adjusted 
to also reproduce the first two orders c\ \ and c-l\- The motivation for also making use 
of the fifth datum, namely the estimate for the coefficient 05,1 = 283, as well as including 
the polynomial terms, is that we wish to avoid fixing the residue of a renormalon term 
by using orders as low as n = 2. However, as we discuss below, our result is surprisingly 
independent of this extra assumption. 

Following the outlined procedure, the parameters of the model ([O) are found to be: 



d^ Y = -1.56- 1(T 2 , 4 R = 3.16, 4 R = -13.5, (6.2) 



10 The additional "1" arises from the anomalous dimension term 7q in eq. (5.11) related to a four- 
quark operator fjsfj . 

11 A similar approach to the heavy quark mass has been put forward in refs. 



C6,l 



C8,l 



C9,l 



ClO.l 



Cl2,l 



3275 1.88 -10 4 3.88 • 10 5 9.19 • 10 5 8.37 • 10 7 -5.19 • 10 8 3.38-10 
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Table 2: Predictions for the Adler-function coefficients cq.i to ci2,i from our model (3.1) for 
B[R](u), employing the estimate for the coefficient c 5 ,i — 283 of eq. ( 3.1C ). 
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Figure 6: Results for D(M^) (full circles) at a s (M T ) = 0.34, employing the higher-order coefficients 
c nj i of table |^ obtained from our model ( |6.l[ ), as a function of the order n up to which the terms in 
the perturbative series have been summed. The straight line represents the result for the Borel sum 
of the series, and the shaded band provides an error estimate inferred from the complex ambiguity. 



d?° = 0.781 , 4° = 7.66 • 1(T 3 . 

The fact that the parameter d^ turns out to be so small implies that the coefficient 
C2 t i is already reasonably well described by the renormalon pole contribution, although it 
was not used to fix the residua. Another implication of this observation will be discussed 
below. With the Borel transform thus determined, we are now in the position to investigate 
the higher-order coefficients. The perturbative series of the Adler function up to ci2,i is 
presented in table ||. We observe that the first negative coefficient is found at the 11th order, 
after which the series retains its sign-alternating behaviour. The successive approximations 
to the (reduced) Adler function D(M^) defined in ( |4.4j ) are displayed in figure |(| which 
shows that the model is well-behaved: the series goes through a number of small terms at 
order n = 4 to 7, the minimal term being reached at n = 5, such that the truncated series 
agrees nicely with its Borel sum. The sign-alternating UV renormalon divergence takes 
over around n = 10. 

To gain further insight into the contribution of a certain renormalon singularity to the 
coefficients Cn,i, m table ^ the relative contributions (in %) for the two IR and the UV 
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Table 3: Relative contributions (in %) of the different IR and UV renormalon poles to the Adler- 
function coefficients 03,1 to 012,1. 
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Figure 7: Results for <5pQ (full circles) and 8 C { (grey circles) at a s (M T ) — 0.34, employing the 
higher-order coefficients c„.i of table || obtained from our model ( |6.l| ), as a function of the order n 
up to which the terms in the perturbative series have been summed. The straight line represents 
the result for the Borel sum of the series, and the shaded band provides an error estimate inferred 
from the complex ambiguity. 

pole are displayed. As can be seen from table ||, already for the third order there is a 
reasonable dominance of the leading IR renormalon at u = 2, and a sizeable but not too 
large contribution of the second one at u = 3. For the next few orders the leading IR pole 
becomes even more dominant as it should. Then, after a region of cancellations between 
the leading IR and UV poles, the leading UV renormalon, which dictates the large order 
behaviour, takes over. Nevertheless, already at the fifth order, the last we have included in 
fixing our parameters, with 16% there is a noticeable contribution from the UV renormalon, 
so that we have sensitivity to this term in the construction of the model. 

Let us now move to the implications of our model (|6.1| ) for the r hadronic width, that is 
for in FOPT and CIPT. A graphical representation of the series behaviours is displayed 
in figure [?]. Like in section ||, the full circles denote our result for 5^q and the grey circles 
the one for 5^1 , as a function of the order n up to which the perturbative series has been 
summed. The straight line corresponds to the principal value Borel sum of the series, and, 
like in figure || for the large-/3o approximation, the shaded band provides an error estimate 
based on the imaginary part divided by n. The order at which the FO and CI series have 
their smallest terms is indicated by the grey diamonds. For CIPT this happens at the 7th 
and for FOPT at the 8th order. The essential conclusion is that the qualitative behaviour 
of the realistic series is determined by the features that were discussed in section |] in the 
context of the p = 2 single-pole model. The cancellations between the Adler function 
coefficients c Wj i and the contributions from the contour integral, g n , soften the divergence 
of the R T series relative to the Adler function series shown in figure |6|, allowing the FO 
series to approach the "true" result around its minimal term, after which the large-order 
asymptotic behaviour takes over. On the contrary, CIPT misses this cancellation and 
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always stays much below the "true" result. If CIPT were used to determine a s from the 
measured value of R T , i.e. 5^°\ a too large value of a s would be extracted to compensate 
for the deficit. We emphasise that the clear preference for FOPT holds even if we discard 
all higher-order terms and use only the series up to n = 4, as can be seen from figure [7|, 
provided the gross features of our ansatz are correct. 12 

Numerically, our central result for the Borel sum is found to be 

6<£l = 0.2371 ±0.0060 i. (6.3) 

Due to the suppression of the leading IR renormalon divergence by the contour integration, 
only one third of the imaginary part arises from the first IR pole, and two thirds from the 
second. This is consistent with the power-suppressed terms in the OPE, where for R T the 
1/M® terms dominate over the terms for the same reasons. The ambiguity that our 

procedure assigns to the "true" value is also of order of the power corrections as discussed 
in section At the order of its minimal term, the FOPT result reads 6^q = 0.2353 with 
a minimal term of 0.0011. Thus, the difference of the Borel sum and the supposedly best 
result for 5pQ is about 1.6 times the minimal term, and the imaginary ambiguity divided by 
7r precisely is able to account for this difference. Thus, our ansatz passes the requirements 
for a sensible asymptotic series and OPE expansion. In the remainder of this section, we 
shall make further checks as to the robustness of the prediction for 5^ from our model for 
B [D] (u) , and we estimate the uncertainty of our value for <5g°g given in fl6.3|) . 

Our first test consists in dropping C5 ; i as an input, and instead using C2,i to determine 
the residua of the renormalon poles. In this case we set since we have only 

four data to fix the parameters. The visual appearance of the corresponding plot for 
<5(°) is indistinguishable from figure [5], so we do not show it again. Also numerically, 
with <5gg = 0.2370 ± 0.0061 i, we have practically the same result as in our main model. 
However, now we are in a position to predict the coefficient C5 1, with the surprising result 



05,1 = 280, almost the same as our naive estimate (3.10) that was used as an input before. 
This outcome was already anticipated by the fact that for our main model ( |6.1| ) , d^ (which 
we now set to 0) turned out to be very small. However, when one breaks down 02,1 into 
the contributions from the different renormalon terms, one finds that the IR pole at u = 3 
is twice as important as the leading IR pole, so the agreement with the previous ansatz 
might be fortuitous. For this reason, we prefer to base our discussion on an estimate of 
C5 ; i, despite the obvious drawback of having to rely on an additional assumption. 

Another modification to corroborate our assumptions and the findings for the main 
model is to add a third IR renormalon pole at u = 4, again only adding one constant rfj" , 
and fitting the parameters to the first five perturbative coefficients. Hence, like in the last 
model, the coefficient C2,i influences the residua of the renormalon poles. The numerical 
outcome is <5gg = 0.2377 ± 0.0065 i, once more only a slight shift compared to (|6.3|). The 
small change in the complex ambiguity only originates from the first two IR poles with 
almost no contribution from the third. Also graphically, no difference would be visible in 
comparison to figure |7|. When inspecting the separate pole contributions, the third IR pole 



12 Further insights into the origin of the difference between <$™ and S^I can be found in appendix B. 
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mainly contributes to the first and second coefficient. For c^ t \, we are left with a small 5% 
contribution, and beyond the third order, the IR pole at u = 4 is completely negligible. 
Thus, our assumptions on the relevant singularities appear to be justified. 

Two further sources of uncertainty in our result fl6,3|) arise from the dependence on 
the estimate of the coefficient 05,1, which is not exactly known, and from neglecting higher 
order running effects beyond (3^, where the /3-function coefficients are not known. In view 
of the fact that 05,1 turned out extremely stable when it was predicted above from our 
model, to estimate the corresponding uncertainty, we think it is a save choice to vary 05^ 
by ± 50%. This results in variations of (± 0.0042, ± 0.0030) in the real and imaginary part 
of <5gg, respectively, with no qualitative change in the behaviours of the FOPT and CIPT 
series relative to figure [7|. 13 To get an idea about the importance of higher order running 
effects, we can set /?4, 62 and 62 to zero, and check how much our result changes. This 
results in a (—0.0021, —0.0011) variation of <5gg. Adding both variations in quadrature, 
we arrive at our final result for <5gs including uncertainties, 

4s = (0.2371 ± 0.0047) ± (0.0060 ± 0.0032) i , (6.4) 

which will be used in the next section for our determination of a s from the hadronic r 
decay width. This should be compared to the perturbative corrections from eqs. ( jOj ) and 
(U), which give 0.2200 (0.2288) in FOPT and 0.1984 (0.2021) in CIPT when the series is 
truncated at n = 4 (n = 5). 

Note that there is no uncertainty due to scale- or scheme-dependence in the usual 
sense, since we are considering a series to all orders that is formally scale- and scheme- 
independent. A question that can be asked, however, is whether our ansatz for the Borel 
transform would still be meaningful and lead to the same results regarding the validity 
of FOPT and CIPT, if the input data c n i, n < 5, were given in another scheme for 
a s than the MS scheme or at another scale. This is not obvious, but it is not clear 
what conclusions should be drawn from this. An arbitrary scheme change can produce 
arbitrary irregularities in the input data, making any attempt to merge low with high orders 
meaningless. Similarly, a scheme change defined by the relation l/a s = l/a^ s + (3\C /(2tt) 
would produce an additional factor exp(— Cu) in the new Borel transform, which changes 
the renormalon residues and therefore the balance between contributions from the different 
leading singularities. With no guidance at hand, we adopt the point of view that the MS 
scheme has proven useful and stable in so many applications of perturbative QCD that 
there is little motivation to consider significant departures. Nevertheless, to acquire an 
estimate of the model dependence of our approach, in the next section we also fit our 
model parameters to the Adler function coefficients which correspond to the expansion 
in the strong coupling evaluated at a scale a s (£M T ), and determine the impact of this 
variation on the a s determination. With such variations, the strength of the leading UV 
renormalon is modified as compared to the lowest IR renormalon pole. For £ < 1 the 

13 This remains true for larger variations of 05,1. However, for 05,1 = or negative, large cancellations be- 
tween the two IR renormalon poles are required in our model to produce such small 05,1 and the assumption 
to fit the low orders with only a few renormalon singularities becomes questionable. 
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leading UV renormalon becomes stronger, leading to an earlier onset of the asymptotic 
regime, while for £ > 1 it is suppressed, which imposes limitations on sensible values of £. 

7. Determination of ct s 

The starting point for a determination of a s from hadronic r decays is eq. ( |1.2|) for the decay 
rate of the r lepton into light u and d quarks. The general strategy for our extraction of a s 
below will be as follows. We concentrate on the sum of the vector and axialvector channels, 
R T y + A, as in this case some of the higher-dimensional operator contributions cancel. A 
more elaborate analysis of the separate channels along the lines of refs. (?], ||, || employing 
moments of the r decay spectral function is left for the future. We expect significantly 
larger perturbative corrections and ambiguities for the moments, along with the enhanced 
condensate contributions, which makes the a s analysis more complicated. We base our 
analysis on FOPT, and FOPT together with the ansatz for higher-order terms of section || 
but do not use CIPT, despite the fact that the CI perturbation series appears to be better 
behaved in low orders. The reason is that the analysis in the previous sections shows that 
the better stability and smaller scale dependence of CIPT reflects the plateau of the CIPT 
curve in figure [?], but the plateau value of the perturbative correction is far from the true 
value. In such a situation, theoretical error estimates based on scale dependence provide a 
very large underestimate of the true uncertainty. 

The first step of the a s analysis consists in estimating the values of the power correc- 
tions 5^P V+A in eq. ( |1.2| ), which arise from higher-dimensional operators in the framework 
of the OPE. The power corrections will be calculated in FOPT. 14 Given these estimates and 
experimental data, we calculate a phenomenological value of 5^ using eq. ( |l.2[) . We then 
determine the value of a s (M T ) by requiring that the theoretical value 6^ eo in FOPT, and 
in our model of section ^ matches the phenomenological value S^ en . Errors are estimated 
by varying all parameters within their uncertainties. 

7.1 Estimate of power corrections 

Before going through the estimation of the power corrections, let us remark that these will 
only be considered for the transversal part of the correlator Ii^ A '^ 1+0 \s). The longitudi- 
nal contribution cW,s+p to R T y + A, arising from U^ A '^ \s) which is related to scalar and 
pseudoscalar correlation functions, will be included later according to a phenomenologi- 
cal approach which was already employed in the determination of \V US \ from R T E7| , [l8| . 
This approach has much smaller uncertainties than making use of the corresponding QCD 
expressions for the scalar/pseudoscalar light-quark correlators since their perturbative ex- 
pansions converge substantially slower than those for the vector /axial vector correlation 
functions. With foresight, we shall use a s (M T ) = 0.3156 ± 0.006 in the numerical estimate 
of the power correction below. 

14 In principle, for quantities other than the purely perturbative corrections to the Adler function, the 
preference of FOPT over CIPT has to be investigated anew. However, the power corrections are already a 
small correction, and for our application the difference in treating the perturbative expansion of the Wilson 
coefficients is irrelevant. 



- 26 - 



The lowest-dimensional power corrections to R T y+A are of dimension-2 and only arise 
from terms proportional to quark masses squared. Detailed expressions for these contri- 
butions were given in ref. (2^, |4t| up to order a 2 , and the next corrections at O(a^) were 
first presented in ref. []50]]. From these expressions we calculate the corresponding 5^ 2 y+A 
in FOPT. Contributions to this quantity arise from two sources: one is proportional to 
(m 2 + to 2 ,) and the other to to 2 . The second is due to internal strange-quark loops and 
only starts at 0(a 2 ). While the perturbative series for the (to 2 + TO 2 ?) term displays a rea- 
sonable convergence, the to 2 one is very badly behaved, such that in total the O(a^) term 
is about a factor of three times the 0(a 2 ) term. To estimate the uncertainty in 5^ 2 v+A^ 
we average the results either including or omitting the O(a^) term, and take the spread as 
the error. For the quark masses we use m u {M T ) = 2.8±0.5MeV, m,d(M T ) = 5.0±0.6MeV 
and m s {M T ) = 97 db 9MeV, which derive from ref. ]51j| . Also quadratically including the 
parametric uncertainties, which albeit play a minor role, yields 



(2) 

m 2 ,V+A 



6' 2 



(3.1 ±8.6) • 10" 



(7.1) 



Although this contribution has a very large uncertainty, it will turn out to be immaterial 
for the a s determination. 

Dimension-4 contributions arise from three possible sources: the gluon condensate 



(aG ), the quark condensate (qq), and mi corrections. (See e.g. ref. [22| and references 



therein for explicit expressions.) The quartic mass corrections are tiny, and thus we shall 
drop them. Furthermore, being suppressed by 1/s 2 because of the weight function in R T , 
the contour integral is only non-vanishing if there are additional logarithms ln(— s), which 
first appear at order a 2 . 15 The explicit expressions for the two contributions including the 
known terms are then found to be 



\g 2 ),v+a 



? (4) 

\qq),V+A 



IItH 



54 if 



a(M, 



2\2 



Mi 



(m u + m d )(qq) 



a(M 2 ) 2 + 



517 
36 



C3 + 



( 3 )kR ( 



(7.2) 



(7.3) 



where k = (ss)/(qq), R s = 2m s /(m u + to^), and we have assumed isospin symmetry 
for the up and down quark condensates. Historically, the standard value for the gluon 
condensate is (aG 2 ) = 0.012 GeV 4 j53], and not much progress has been made since then. 
Therefore, it will be employed as our central value. As we only know the leading term for 



, to be conservative, we assign a 100% uncertainty. For S^L^y +A , the required 



x(4) 

°{G 2 ),V+A 

quark condensate can be calculated from the GMOR relation [53, 54], with the result 
(qq)(M T ) = — (272±15MeV) 3 . With 40%, the next-to-leading order correction is large, 
but still has perturbative character. Thus, we include this term and take its size as an 



estimate for the missing higher orders. Further employing k = 0.8 ± 0.3 [54], and the light 
quark masses from above, we obtain 



r(4) 

\g 2 ),v+a 



(3.3 ±3.3) • 10" 



\qq)y+A 



(-4.9 ±6.2) • 10" 



(7.4) 



J This is yet another manifestation of the suppression of the leading IR renormalon at it = 2. 
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It may be remarked that the upper range of S^ 2 \ V+A just corresponds to the complex 
ambiguity of the u = 2 renormalon pole divided by tt in our all-order ansatz. It is gratifying 
to observe that they are of the same order of magnitude, giving support to the size of 
this term. Nevertheless, also the dimension-4 contributions only play a minor role in the 
determination of a s . 

At dimension six, S^y, A receives contributions from the three-gluon condensate 
(g 3 G 3 ), 4-quark operators, and lower-dimensional operators times appropriate powers of 
quark masses. The coefficient function of (g 3 G 3 ) vanishes at the leading order and is there- 
fore suppressed. Also the operators which get multiplied by quark masses only give very 
small contributions. Thus we neglect these two types of contributions and concentrate on 
the 4-quark condensates. At leading order, one is confronted with three four-quark opera- 



tors, while at higher orders many more are generated y, 55, 56 1. As it appears impossible to 
determine all required condensates from phenomenology, the so-called vacuum-saturation 
approximation (VSA) had been proposed p^j. This reduces all 4-quark condensates to 
squares of the quark condensate. With this simplifying assumption, the leading dimension- 
6 contribution takes the form 

<8L.V*. = - f A « ^ = (" « ± 2.9) ■ ID- , (7 . 5) 

where we assume p = 2 ± 1 for the numerical estimate. Our reasoning for this is as follows: 



VSA is incompatible with the scale and scheme dependence of the 4-quark operators [56 



57]. Therefore, it does not make sense to include the next-to-leading order corrections. 
Conventionally, one introduces the parameter p, which comprises the violation of the VSA. 
A rough idea about the size of the parameter p can be gleaned from phenomenological fits 



of the dimension-6 contributions to IT^ d , H"^, II n( ^~ , as well as H ud |7|, ||, ||, [58| (and 
references therein). Comparing the phenomenological fits with the VSA, values for p from 
one to about three are found, which motivates the chosen range. At any rate, "/J q ^ y+A 
constitutes the dominant power correction. 

We also include a crude estimate for the still higher-dimensional operator corrections. 
Not much is known about these contributions, apart from the fact that they should be 
suppressed as compared to "^Lgg) y+A' smce they carry at least two more powers of 1/M T . 
Inspecting the fits of refs. j?], ^ || where also a dimension-8 contribution was included, one 
infers that roughly it could be of order 10~ 3 . Therefore, we have added 

Cv + a = (0 ± 1) ■ ID" 3 (7.6) 

in our total estimate of the uncertainty for power corrections. There are potentially also 
(short-distance) instanton contributions to the r hadronic width. The leading contribution 
of this type is a rapidly increasing and uncertain function of a s (M T ). For a s (M T ) = 0.32 it 
has been estimated to contribute 2 • 10~ 3 to S^y +A Nonetheless, we do not include 
a further power correction uncertainty for this term. 

To complete our summary of power corrections to R T , we still have to compute the 
longitudinal contributions which arise from scalar and pseudoscalar correlators. Because 
the perturbative series for these correlators do not converge very well, here we shall follow 
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the approach of refs. [47, [Ig]. The contribution from the scalar correlator is suppressed by 
a factor (m u — md) 2 , and it can be altogether neglected. The main idea then is to replace 
the QCD expressions for the pseudoscalar correlator by a phenomenological representation. 
The dominant contribution to the pseudoscalar spectral function stems from the well known 
pion pole, giving 

flMl[ MgV 
M T 4 1 A/2 ) 



= ~ ^ ^ I 1 - TH I > ( 7 - 7 ) 



plus small corrections from higher-excited pionic resonances. Repeating the analysis of 
section 3 of ref. E?]] and updating the input parameters, we find 



<W+P = (- 2.64 ± 0.05) • 10~ d . (7.8) 

The uncertainty in (|7.8| ) has been estimated by setting the higher-resonance contribution 
to zero, which clearly demonstrates that the pion pole is dominant. Collecting all contri- 
butions, and adding the errors in quadrature, we arrive at our total estimate of all power 
corrections: 

5 PC = (-7.1±3.1)-10~ 3 . (7.9) 

Our value (^^) is consistent with the most recent fit to the r spectral functions performed 
in ref. ||. 

As a matter of principle, the OPE of the correlation functions in the complex s plane, 
even when integrated over a suitable energy interval could be inflicted with so-called 
"duality violations" [6C[. In our case, these would come from the contour integral close 



to the physical region, and even though suppressed by a double zero, could lead to addi- 
tional contributions or uncertainties. Within a model, these contributions were recently 
investigated for hadronic r decays and in ref. j|] the model was fitted to the r spec- 
tral functions, with the finding that possible additional contributions are below the 10 -3 
level. In view of these results, we shall omit possible duality violating terms, but in future 
analyses which perform a simultaneous fit of higher order OPE contributions, in analogy 
to |J, it might be worthwhile to include them. 

7.2 a s analysis 

Employing the value R T y+A = 3.479 ± 0.011, which results from eq. (|1.1| ) in conjunction 



with R TiS = 0.1615 ± 0.0040 @, as well as \V ud \ = 0.97418 ± 0.00026 §|, from eq. fllj) 
the phenomenological value for <j(°J can be derived after accounting for the electroweak 
corrections and subtracting the power correction given in eq. (|7.9|): 



5^ en = 0.2042 ± 0.0038 exp ± 0.0033 PC = 0.2042 ± 0.0050 . (7.10) 

By far the dominant experimental uncertainty is due to Rry+A- In the second error, we 
have also included the ones from 5"ew an d S' EW which should be considered theoretical 
(though they are not power corrections). The final step in the extraction of a s (M T ) now 
consists in finding the values of a s for which the phenomenological value <^p^ en matches the 
theoretical prediction, either from the finite order series, or from our model to all orders. 
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Let us begin by analysing the impact of using FOPT when the perturbative expansion 
is employed up to and including the fifth order. Besides the exactly known coefficients, we 
make use of our value ( 3.1C| ), 05,1 = 283, and to estimate the corresponding uncertainty, we 



have either removed or doubled the fifth order term. A further theoretical error is added 
to account for the residual renormalisation scale dependence of the fixed order series. We 
estimate this by expressing the FO series in terms of a{fx 2 ) rather than a(M^) and by 
varying /i between 1 GeV and 2.5 GeV. Going through this procedure, and keeping the 
errors separate to clearly see the importance of the various contributing uncertainties, we 
obtain: 

a s {M T ) = 0.3203 ± 0.0032 exp ± 0.0028 PC ± 0.0027 C51 1 HoeI (scale) 

= 0.320 ± 0.003 CX p toUl (th) = 0.3201 • ( 7 - n ) 



Comparing our result ( [7.11 ) with other determinations of a s (M r ) from hadronic r decays 



H [l8| that include the exact fourth order term, we observe that our value is smaller by 
0.012 [18] and 0.024 Q, up to twice the estimated error. This should not come as a surprise 
since the analysis of ref. || relied on the use of CIPT, which increases a s (M T ), see figure [7|, 
while in ref. [18] an average of FOPT and CIPT is performed. The error in eq. ( 7.11| ) 



is dominated by the scale error, which is significantly larger than in CIPT. 16 However, 
as argued above, we find that the small scale-dependence in CIPT is misleading when 
interpreted as a measure of the theoretical uncertainty. 

We evolve a s (M T ) to the Z-boson mass scale, using the /3-function known to four 
loops 63, 64], and the matching coefficients at flavour thresholds up to order [65]. 17 



The central values of the flavour thresholds /i* = 3.729 GeV and \x* h = 10.558 GeV are 
taken to correspond to the physical thresholds of open D and B meson production. In 
the matching coefficients also the charm and bottom quark masses are required, which 
we assume to be m c (m c ) = 1.28 ± 0.05 GeV and m b {m b ) = 4.20 ± 0.05 GeV @ in the 
MS-scheme. The uncertainties resulting from the evolution are estimated as follows: like 
in section ^, we compare to consistent 3-loop running, i.e. setting (3^ = and removing 
the O(a^) coefficient in the matching relation, and we vary the quark masses as well as 
the flavour thresholds in the ranges \x* c j1 < /i c < 2/i* and nl/2 < ^ < 2/i£. Adding all 
uncertainties in quadrature, we arrive at: 

a s (M z ) = 0.1185 ± 0.0004 cxp jl°;°°ig (th) ± 0.0002 evol 

= 0.1185±g : ggS. (7-12) 

Again, this value is lower than previous results from hadronic r decays, but it is in perfect 
agreement with the global averages of a s (Mz) |l], |2|. 



16 The scale variation has a minimum at the scale [i ~ 1.22 GeV, which might therefore be considered as 
an "optimal" scale. It is amusing to note that evaluating a s at this scale, we obtain a s (M T ) = 0.3151, very 



close to the improved result (7.13). 

17 The evolution has been performed independently with the Mathematica package RunDec and a 
private routine coded by one of us (MJ), finding complete agreement. 
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To further improve our determination of a s , we now include the higher order pertur- 
bative terms according to our model of section 0. To this end, we need to find the value 



of a s such the value of the Borel sum 5g°g matches the phenomenological value ( T.1C| ). 



is 



Again keeping the contributing uncertainties separate, we find: 

a s {M T ) = 0.3156 ± 0.0030 cxp ± 0.0026 PC ± 0.0025 cg A ± 0.0011^=0 ± olooio ( scale ) 

= 0.3156 ± 0.0030 cxp ± 0.0051 th = 0.3156 ± 0.0059 . (7.13) 



Besides the uncertainties already discussed previously, in eq. ( 7.13j) we have also included 



a scale/model uncertainty according to the reasoning put forward at the end of section |6|. 
As the scale variation, we chose 0.5 < £ 2 < 1.5, where fi = £M T . (The corresponding 
scale interval is 1.26 GeV < \i < 2.18 GeV.) The employed range in £ is dictated by the 
observation that for lower scales the contribution of the leading UV renormalon is enhanced 
as compared to the leading IR renormalon pole, while for larger scales it is suppressed. At 
the lower end for £, the contribution of the leading UV renormalon to the coefficient C3 i is of 
a similar size than the leading IR renormalon, while at the upper end, we are only left with 
a mere 3 % contribution of the first UV renormalon to the coefficient 05,1, practically losing 
the sensitivity to this contribution. Therefore, beyond the used range for £, our model 
ceases to make sense. For given £, we determine a s (£M T ) such that the value (|7.1C| ) is 
obtained and then evolve back to the scale M T . In quoting the final uncertainty of a s (M T ), 



in ( 7.13 ) we have used the larger scale variation in order to have a symmetric final error. 
The dependence of the successive perturbative approximations to R T on the choice of £ is 
shown in figure ||, upper panel. The lower panel shows the corresponding approximations 
in CIPT, which are again seen to lie below the "true" result for any reasonable value of £. 
Evolving the result ( [7.13 ) to the Z-boson mass scale, we arrive at our final value for 



a s (M z ): 

a s (M z ) = 0.11795 ± 0.00038 cxp ± 0.00063 th ± 0.00020 evo i 

= 0.11795 ±0.00076. (7.14) 



This result is slightly smaller than eq. ( 7.12| ) based on the FOPT up to the fifth order as 



could be anticipated from figure 0, where one observes that 5^ at 0{a\) is smaller than 



the full Borel sum. We consider eq. ( 7.14 ) as our best estimate of the strong coupling in 
the MS scheme from hadronic r decays. 



8. Conclusions 

Hadronic t decays provide an especially clean environment for the study of QCD effects 
and in particular the determination of QCD parameters. Of prime interest in this respect 
is the QCD coupling a s . Still, due to the relatively low scale M T , an adequate control 
over the perturbative series should be achieved for such applications. Besides explicitly 

We do not include a separate error from the ambiguity of 8$ (at a 3 (M T ) = 0.3156: 4s = 0.2042 ± 
0.0029 i), since it is subsumed in the error of <5pc, which is several times larger. 
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Figure 8: Results for S^q (full circles, diamonds, triangles) and 5^ (grey circles, diamonds, trian- 
gles) for the model of eq. ([O]) as a function of the order n up to which the terms in the perturbative 
series have been summed. The parameters of the model are fit to the first five coefficients of the 
Adler function expanded in a s (£\/s) as described in section [| The value of a s (£M T ) is determined 
by requiring that the Borel sum equals 6^ en = 0.2042 in each case. 



computing terms in the perturbative expansion, progress in this direction can be attained 
by inspecting the renormalon divergence structure that the perturbative series should have 
on general grounds. 

A long-standing question in the interpretation of the QCD correction to R T is a numer- 
ical discrepancy between two ways of performing the renormalisation group improvement, 
namely CIPT and FOPT. While CIPT resums running effects of a s along the complex 
integration contour which are known to be large, FOPT performs a consistent expansion 
in a s at each loop order. The CIPT sum is significantly below the FOPT sum requiring 
a larger strong coupling to reproduce the experimental R T . Resolving this discrepancy 
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has become a major issue for improving the accuracy of the a s determination from R T , 
in particular since recently the fourth order coefficient of the series expansion has been 
computed [O], and further terms cannot be expected to be available any time soon. 

CIPT has long been considered as the method of choice due to its apparent better con- 
vergence and smaller renormalisation scale dependence of the truncated series as compared 
to FOPT. However, the argument that FOPT should be discarded, since the expansion of 
the running coupling along the contour produces a series of finite radius of convergence that 
is avoided in CIPT, does not take into account that both the CIPT and FOPT series have 
zero radius of convergence anyway due to factorial divergence in higher orders. This leads 
to interesting cancellations between the Adler function coefficients and running coupling 



effects |pL7 |. Having a sufficiently large number of exact Adler function coefficients at our 
disposal, and without prospects of further improvements by exact computations, makes it 
timely and possible to attempt to merge the exactly known terms with what is known on 
general grounds about large-order behaviour, and to reinvestigate the conceptual issues of 
the FOPT/CIPT comparison. 

We investigated several toy models to find out under which circumstances FOPT or 
CIPT provide a better approximation to the "true" result, which we may define as the 
Borel sum of the series, since power corrections to R T are small. Two extreme cases can 
be singled out. If the Adler function coefficients are set to zero beyond a certain order, 
CIPT becomes exact, while FOPT performs large oscillations around the exact result [19]. 
On the other hand, if the Adler function series is dominated by an infrared renormalon 
at u = 2, large cancellations occur in the R T series, which are only manifest in FOPT. In 
this case, FOPT accounts well for the Borel sum, while CIPT stays systematically below, 
despite an apparent better convergence and smaller scale dependence! 

Taking these lessons we constructed a realistic ansatz for the entire Borel transform 
of the Adler function such that the exactly known coefficients plus an estimate of the fifth 
order coefficient are exactly reproduced, while including knowledge about the leading three 
renormalon singularities. Given the regularity of the exactly known low-order terms, the 
weight of the three contributions was fit to the four- to six-loop terms. With our central 
model of eq. ([0]) we were able to achieve a very coherent picture of the known and higher 
order coefficients. Our main conclusion from a study of the realistic model is that, given the 
particular features of the R T series, FOPT provides the better approximation, in general 
and at fourth order, than CIPT and is to be preferred on these grounds. In this respect 
the full QCD case resembles the large-/?o approximation discussed in ref. p7f| . 

We believe this conclusion to be valid in QCD beyond particular models, since it is 
based on generic properties of the R T perturbation series - although it must be admitted 
that statements about the accuracy of FOPT and CIPT, and the validity of perturbation 
theory in higher orders in general, are always open to some amount of speculation. The 
main point is that conventional wisdom that favours CIPT is based on the first three orders 
in perturbation theory when the running coupling effects in g n dominate over the Adler 
function coefficients c n ^. This dominance weakens as n grows, as can be seen from the 
known exact coefficients, and the situation is expected to reverse beyond n = 5. The 
fact that the c n> i ultimately diverge, cannot be ignored. This as well as the crucial 1/ra 2 
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cancellation in the sum c n> i + g n for the leading non-sign-alternating component of the R T 
perturbation series are model- independent consequences of QCD. 

Making use of these results, two determinations of the strong coupling a s from R T y + A 
were presented in section |7[ First, we performed a conventional FOPT analysis based 
on the fourth-order result plus an estimate of the fifth order, whose size is taken as an 
uncertainty. The resulting value for a s (M T ) can be found in eq. ( [7,11| ). Evolving to Mz, 
we obtain 

a.(M z ) = 0.1185t°;°oo9 (FOPT), (8.1) 



where the error is dominated by residual renormalisation scale dependence; see eq. ( 7.12| ) 



This value is lower than those presented in refs. || 15] mainly because we propose to not use 
CIPT as the result of our study. Our preferred and best result is obtained by incorporating 
available structural information on the large-order behaviour as embodied in the model of 
eq. (|6.l| ). The resulting value for a s (M T ) can be found in eq. (|7l3b - Evolving it to M z , 
the main result of our paper reads: 

a s (M z ) = 0.11795 ± 0.00038 exp ± 0.00063 th ± 0.00020 cvol 

= 0.11795 ±0.00076. (8.2) 

Our approach towards the investigation of the higher-order behaviour of perturbative 
series employed in this work can certainly also be applied to other quantities of interest, 
like the scalar correlation function, mass squared corrections to R T , or moments of the 
spectral functions. It should be interesting to see what can be said about the issue of 
renormalisation group resummation when contour integrations are involved in these cases. 
We shall return to this question in the near future. Finally, our results can be useful in 
other places where the QCD Adler function plays a role. 
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Appendix A: Borel integral for renormalon poles 

In this appendix, we provide a few useful expressions for the Borel integral of generic 
renormalon poles, and the corresponding complex ambiguity in the case of IR renormalons. 
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Let us begin with a general UV-renormalon pole which we assume to have the form: 



B[R^](u) 



(A.l) 



+ up ' 

where u = Pot with Po = Pi/(2tt), p £ N and 7 £ R + . Then, the corresponding Borel 
integral is given by 



R^ v {a) = d 



uv 



-t/a 







+ Potp 



dt, 



(A.2) 



which after the substitution t = az — p/ (3q can be expressed in terms of the incomplete 
T-function: 

00 



R% V {a) = d 



•uv 



a 



7 



. gP/(/3oo!) 



z 'e 



dz = d. 



uv 



a 



p/(A>a) 



eP/0fc«)r(l- 7 ,^). (A.3) 



Expanding the UV-renormalon pole ansatz (A.l), and performing the Borel integration 
term by term, yields the corresponding perturbative expansion: 



o 



ra+l 



(A.4) 



For the general IR-renormalon pole, let us assume the generic form: 



4 R 



p JW ~ (p-up ' 

where again u = Pot, p£N and 7 G M + . The corresponding Borel integral is given by 



(A.5) 



i?i R (a) = d 



■1R 



-t/a 



p 7 (p-A)t)T 



dt . 



(A.6) 



and can be expressed in terms of the exponential integral function E 7 (z). However, because 
of the singularity at t = p/Po and the branch cut for t > p/po, we still have to specify the 
defining integration path. We shall integrate on the real axis up to t = p/Po — £, then on 
a semicircle of radius e in a clockwise, or anti-clockwise direction, and finally again along 
the real axis from t = p/Po + £ up to infinity, either above or below the cut, taking the 
limit e — > in the end. Performing all integrals, this leads to the expression: 



r(i- 7 ) . (A.7) 



(_1)±7 _ (_l)sig(Im[a]b 



In the above equation, (— l) 2 should be interpreted as exp(i7r2i), and the function sig(z) 
represents the sign of z, with the additional definition sig(O) = 1. The imaginary ambiguity 
of i?p R (a) can be readily computed from eq. (A.7), and for a > turns out to be: 

d lR 

Im [i?p R (a)] = ± -fy sin(7r 7 ) T(l - 7) a 1 " 7 e - p/(A)a) . 
Po 



(Ai 
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We can also straightforwardly obtain the perturbative expansion of R^(a), which reads: 

n=0 



jIR oo / R \ n 

^ R(Q) = ^? r(n + 7) (f)"" +1 - (A ' 9) 



These results can be used to construct analytic expressions for our models, which can 
all be decomposed into polynomials and a sum of poles of the above form. 

Appendix B: Adler function in the complex s-plane 

Further understanding of the origin for the difference of 5^ as calculated within FOPT or 
CIPT can be gained by inspecting the reduced Adler function D(s) on the circle s = M 2 e 1 ^ 
in the complex s-plane. From the general perturbative expansion ( [2.10 ) it is clear that the 



ambiguity in the choice of the RG resummation already exists in this case. Analogous to 
the discussion of section || for R T , FOPT is defined by employing the constant scale choice 
/x 2 = M 2 , while CIPT corresponds to a resummation of the series on each point of the circle 
separately through the variable choice /i 2 = — s = — M 2 e J</3 . The perturbative expansions 
for D as a function of the angle ip in the two cases take the form: 

oo n 

n=l k=l 

oo 

Dab) = E c «.i <~ M r e ^)" • ( R2 ) 



n=l 



In figure |9], we display a graphical account of our numerical results for the real parts 
Re[DFo(p)] and Re[-Dci(¥>)] m the range < <p < 2tt and for a s (M T ) = 0.3156. We have 
drawn four different lines corresponding to a truncation of the perturbative series at the 4th 
(dotted line), 5th (dashed-double-dotted line), 6th (dashed-dotted line) and 7th (dashed 
line) order. In addition, the solid line corresponds to a resummation of the perturbative 
series according to the Borel sum of the physical model for the Adler function introduced 
in section ||. The required perturbative coefficients cq : i and cj : i for this model are given in 
table |2| The reason for showing the summation of the series up to n = 7 lies in the fact 
that the smallest summand of both series is found in this range. 

Various observations can be made on the basis of figure ||. First of all, as can also be 
seen from eqs. ( |B.1| ) and ( [B.2D , at the euclidian point p = n on the negative real s-axis, 
Dpo an d Dqi are identical. The closest approach of the partial sums to the Borel sum for 
<p = 7r is reached at the 6th order (compare figure ||) . Close to the minkowskian region (the 
positive real s-axis), that is ip near zero or 2tt, FOPT converges rather slowly, and only 
at the 7th order, the series approaches the shape displayed by the Borel sum. The slow 
convergence of FOPT close to the minkowskian region has already been noted in ref. ||, 
and was taken as an argument that CIPT should be preferable. Indeed, as is obvious from 
the lower plot of figure ||, CIPT converges very nicely for all angles and mostly even better 
for ip not near tt than at the euclidian point. However, the resulting sums do not at all 
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Figure 9: The real parts Re[Z?Fo(<p)] and Re[Dci(tp)] °f ec t s - ( B.l ) and (|B.2| ) for a summation of 
the perturbative series up to four different orders, namely the 4th (dotted line), 5th (dashed-double- 
dotted line), 6th (dashed-dotted line) and 7th (dashed line) order. The input for the QCD coupling 
was chosen to be a s (M T ) = 0.3156. Also displayed as the solid line is the Borel sum Re[Z?Bs( ( / ? )] 
according to our physical model presented in section |6[ 



approach the shape of the Borel sum. Thus, except at a few isolated points, CIPT for 
the Adler function in the complex plane exhibits a similar problem as for R T , shown in 
figure |7], that the series exhibits apparent convergence and stability, but the resultant value 
is much further from the "true result" than any theoretical error estimate would suggest. 
For the Adler function CIPT lies below the Borel sum in the central region and above 
elsewhere, while for R T the weight function (1 — x) 3 (l + x) with x = e tip in the integration 
over the angle ip implies that 5q^ is below the Borel sum 5g°g as found in section ^. On the 
other hand, FOPT turns out to be much closer to the Borel sum, since the problematic 
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region close to the minkowskian axis is strongly suppressed in R T by the triple zero of the 
kinematical weight function at s = The problems with CIPT for the Adler function 
in the complex plane arise essentially for the same reason as for R T . In FOPT at each 
order we observe sizeable cancellations between the independent coefficients c n .i and the 
running effects which are missed by CIPT, although here, due to the angular dependence, 
the situation is less transparent. 
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